m 
O 
O 
(N 



OS 



Thermodynamical fluctuations in optical mirror coatings 

V. B. Braginsky and S. P. Vyatchanin 
Physics Faculty, Moscow State University, Moscow 119899, Russia 
e-mail: vyat@hbar.phys.msu.su 

February 6, 2008 



Abstract 

<D ; 

^Sermodynamical (TD) fluctuations of temperature 
i^ mirrors may produce surface fluctuations not only 
ttfrough thermal expansion in mirror body pQ but also 
^Jg:ough thermal expansion in mirror coating. We ana- 
Jijde the last "surface" effect which can be larger than 

t first " volume" one due to larger thermal expansion 
fficient of coating material and smaller effective vol- 
,^llne. In particular, these fluctuations may be impor- 
1j*Mit in laser interferometric gravitational antennae. 

o 

4- 1 Introduction 

^^itstanding experimental achievements in quantum 
pptics and high resolution spectroscopy within recent 
Tfo*r decades are substantially due to creation and de- 
velopment of new technologies. The coating of mirrors, 
[the deposition of many thin dielectrical layers on 
J^e mirror surface is likely to be one of the most impor- 
ftetnt technology in this part of experimental physics, 
^ftgday this technology provides the reflectivity R of 
s^Ccfi coating very close to unity (at the beginning of 
^er era the value of (1 — R) was about several percent 
lile at present (1 — R) ~ 10~ 6 0E]). There is reason 
i hope that in the not too distant future the value 
j(l — R) will be close to 10~ 9 . This improvement 
iJa)particular means that measurement using squeezed 
quantum states will become routine as the squeezed 
ite lifetime is longer if losses are smaller. The rise 
'■^squeezing factor provides the possibility to increase 
ositivity with the same number of "used" photons 

The value of (1 — R) is of great concern in terres- 
trial gravitational wave antennae (projects LIGO and 
VIRGO, e.g. see 00 IE])- These antennae can be re- 
garded as extremely sensitive spectrometers which goal 
is to detect very small deviations Acufp/lufp of one 
Fabry-Perot optical resonator mode eigen frequency 
in comparison with another one. The amplitude of 
these frequency relative deviations which differs from 
the metric perturbations amplitude only by the factor 
of two, is quite small. The planned resolution for the 
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Figure 1: Illustration for qualitative consideration 
of temperature thermodynamic fluctuations producing 
the thermoelastic noise 
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and for the stage 



stage LIGO-I is Acufp/cufp 
LIGO-II is one order smaller. 

The above considerations show that the value of 
(1 — R) decrease is important for further development 
of many optical measuring devices types. However, 
to our knowledge there is no deep analysis of possible 
additional noises which can appear specifically due to 
coating. In this article we present the analysis of a 
noise source in the coated mirror, which deserves seri- 
ous attention. Thermodynamical (TD) fluctuations of 
temperature in mirror coating can produce additional 
fluctuations of surface through thermal expansion in 
mirror coating. This noise may be larger than analo- 
gous one in mirror body due to larger coating materials 
thermal expansion coefficient and smaller effective vol- 
ume. 

2 Semi-qualitative consideration 

"Bulk" TD fluctuations. We have shown in £Q 
that thermodynamical fluctuations of temperature in 
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mirrors are transformed due to thermal expansion co- 
efficient oc.b = (1/l)(dl/dT) into additional noise which 
can be a serious "barrier" limiting sensitivity. This 
noise (it is also called as thermoelastic noise) produces 
random fluctuational "ripples" on surface. This noise 
has a nonlinear origin as nonzero value of ab is due to 
the lattice anharmonicity. 

Our main interest here is fluctuations of averaged 
displacement X(t) of mirror face surface. It is the nor- 
mal displacement v z of the mirror face averaged over 
the beam spot Gaussian power profile: 



x = -L 

7tr; 



o ■> 
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e- r2/r °v z (r, 4>)rd4>dr, (1) 



(to is the beam radius). The spectral density of dis- 
placement X can be presented for half infinite medium 
as follows: 



'bulk ( 



8 ct 2 k B T 2 (1 
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(2) 



Here 1<b is the Bolzmann constant, T is temperature, 
"v is Poisson ratio, k is thermal conductivity and Cv 
is specific heat capacity per unit volume. This result 
has been refined for the case of finite sized mirror by 
Liu and Thorne |12| . however, the variation from our 
result is only several percent for typical mirror sizes, 
and hence we use more compact expression @ here. 

This physical result can be illustrated using the fol- 
lowing qualitative consideration. We consider the sur- 
face fluctuations averaged over the spot with radius To 
which is larger than the characteristic diffusive heat 
transfer length Tt 
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Cvcu ' 
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(3) 



where a> is the characteristic frequency (for LIGO it is 
about w ~ 2n x 100 s _1 ). Since in fused silica Tt — 
3.9 x 10~ 3 cm and in sapphire Tt — 1 .4 x 10~ 2 cm 
the condition Tt <C to ~ 6 cm is fulfilled (material 
parameters are presented in Appendix lB|l . 

We know that in volume ~ r T the length variation 
due to temperature TD fluctuations AT is about 



Axt = abATrT — ocbtiM 
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(4) 



The TD temperature fluctuations in such volumes can 
be considered as independent ones. The number of 
such volumes that contribute to surface fluctuations 



is about N 



5 / r 3 



and hence the displacement X 



averaged over the spot with radius To consists of these 
independently fluctuating volumes displacements sum, 
and is approximately equal to 



x, AXT 



a b r T1 



/k B T 2 
CvTq 1 ' 



(5) 



Comparing JSJ with J2J) we can see that X 



S^ lh Aw (correct to the multiplier of about unity) 
if one assumes Atu ~ cu. This fact confirms that our 
semi-qualitative consideration is correct. 

TD fluctuations in thin layer. One can consider 
the surface fluctuations of half space covered by thin 
layer, produced by TD temperature fluctuations as 
consisting of two parts: "bulk" part depending on mir- 
ror body thermal expansion coefficient <X\, (which has 
been already calculated in JT]) and "surface" part due 
to temperature fluctuations in layer with thickness d 
and effective thermal expansion coefficient a = ot\ — ab 
(here <x\ is thermal expansion coefficient of layer ma- 
terial). 

It is important to note that layer thickness d is 
much smaller than the diffusive heat transfer charac- 
teristic length rj (for optical coating d < 10 urn) and 
therefore the following condition is valid: 



r > r T > d 



(6) 



Therefore we may consider temperature fluctuations in 
our layer to be the same as in layer with thickness Tr> It 
means that TD temperature fluctuations in layer does 
not depend on thickness d. In this case we propose the 
following semi-qualitative consideration: 
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(9) 



Formula JSJ is in good agreement with rigorous formula 
for spectral density Sat(c^) of averaged temperature 
AT in layer obtained in j^j 



S AT (u>) 



V2k B T 2 = V2k B T 2 
nr\\J pCkcu ttCvTqTtuj ' 



(10) 
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and one can estimate the spectral density S 
of surface fluctuations averaged over spot with radius 
To as following 



-TD est l 
layer ^ 



(X 2 & 2 Sat(w) 



v / 2a 2 d 2 k B T 2 
7TpCrQr T cu 



(11) 



Below we will show this formula to give the result cor- 
rect to the multiplier of about unity. 

3 Analytical results 

Using Langevin approach (see also Appendix A) we 
have calculated the field of temperature TD fluctua- 
tions approximating mirror as a half infinite medium 
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with assumption that all material parameters of layer 
and half space are identical with exeption of thermal 
expantion coefficients — (Xb for half space and <X\ for 
layer. These temperature fluctuations produce surface 
deformations of half infinite elastic space covered by 
thin layer. These deformations can be calculated using 
elastic equation (|31[1 (presented in Appendix A). 

We use quasistatic approximation for elastic prob- 
lem because time sound requires to travel across the 
distance To is much smaller than characteristic time 
~ 1/cu. We also assume the layer to be homogeneous. 
We can divide our problem into two problems: 

a) The "bulk" problem. The surface fluctuations 
are produced by TD temperature fluctuations in half 
space covered by layer with different elastic constants. 
The layer and half space have the same thermal ex- 
pansion coefficient at,- It corresponds to the problem 
solved in p. 

b) The " surface" problem. The surface fluctuations 
are also produced by TD temperature fluctuations in 
half space covered by layer. The effective layer thermal 
expansion coefficient oc = oc\ — ctb 7^ 0, and there is no 
thermal expansion in half space. 

The qualitative considerations in previous section 
make clear that TD temperature fluctuations produc- 
ing fluctuations in "surface" and "bulk" problems does 
not virtually correlate and can be calculated indepen- 
dently. 

For "surface" problem we assume that additional 
elasticity produced by layer is much smaller than mir- 
ror body bulk elasticity (due to small thickness of the 
layer) and we take into account only stresses produced 
by layer due to nonuniform temperature distribution 
u(f). Then one can take equation (|31|l and just substi- 
tute into right part the term 



«Vu => adVu6(x — e), 
a = ai — <Xb ■ 



(12) 
(13) 



(Recall that here we assume identity of material con- 
stants for layer and half space.) As a result we obtain 
spectral density of surface fluctuations (see details in 
Appendix |A"|) : 



V2(1 + -v) 2 a 2 d 2 k B T 2 



7T TqVkCvUJ 



(14) 



Here "V is Poisson ration. Note that this formula differs 
from semi-qualitative one (|ll|l only by multiplier of 
about unity. 

The fact that S£° r (tu) > S£°° st -(u>) can be quali- 
tatively explained by the speculation that thermal ex- 
pansion in layer also causes stresses in mirror body un- 
der layer: the material under layer additionally "swells 
up" ("distend"). The similar effect one can observe 
when heated bimetallic plate is bending. 



We have generalized the formula l|14fl for the case 
of different material parameters of layer and substrate 
deriving it from fluctuation-dissipation theorem (FDT) 
(see details in Appendix lB|) : 
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(15) 
(16) 
(17) 



_ Et(1 -2v b 
(1 -v0(1 +^b) E b (1 -"Vi) 

Here Eb, Ei are Young modulus, subscripts b and \ 
refer to parameters of half space and layer correspond- 
ingly. This result was earlier obtained in |1U| . 

We can further generalize the result l|15l - 117|) for 
multilayer coating with alternate material constants 
which differs from material constants of mirror body. 
Let the multilayer coating consists of IM alternating se- 
quences of quarter-wavelength dielectric layers. Ev- 
ery odd layer has the refraction index rii , thickness 
di = A/4rii (A is the incident light wavelength) . Every 
even layer has parameters TI2, d2 = A/4n.2 correspond- 
ingly. Then formulas 1)151 117fl can be used with the to- 
tal coating thickness and effective expantion coefficient 
calculated using the following formulas (see details in 
Appendix [BJ) : 



a 



/VI 

di + d 2 
N(di + d 2 ) 



a b d 2 



d 2 



A 2 



(18) 
(19) 



Here factors A] and A2 are calculated using l|17|) with 
substitution material parameters for odd and even lay- 
ers correspondingly. 

It should be emphasized that negative influence of 
thermoelastic fluctuations in layer in finite sized mir- 
ror can be larger than in half infinite space model due 
to "bimetallic" effect: the temperature fluctuations in 
layer should additionally cause mirror bend through 
thermal expansion. Using FDT approach developed 
by Liu and Thorne 12 we have calculated numerically 
the following coefficient 



Cfsm — 



cTD, finite test mass 
layer 

cTD, infinite test mass ' 
layer 



(20) 



For design planned in LIGO-II the test mass is 
manufactured from fused silica with TQ2O5 + S1.O2 
coating and has following parameters R = 19.4 cm, 
H = 11 .5 cm, To = 6 cm. For this case we estimate 



C* 



1.56 



(see details of rather combersome calculations in Ap- 
pendix EJ). But if R 3> H then the value of Cftm may 
be substantialy larger. 
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4 Numerical estimates 

High reflectivity of mirrors is provided by multilayer 
coating which consists of alternating sequences of 
quarter-wavelength dielectric layers having refraction 
indices n.i and n.2. The frequently used pairs are 
Tq 2 5 (ui ~ 2.1) and Si02 (n 2 ~ 1.45): namely 
these coatings are used in LIGO The total typ- 

ical for high finesse mirror number of layer pairs is 
N ~ 19 (this value we used for estimates below) so 
that total thickness of TQ2O5 + S1.O2 coating is about 
N(di +d 2 )~6x 1CT 4 cm. 

To estimate TD temperature fluctuations in in- 
terferometric coating of mirror we use formula (|15|l 
with substitutions {151 - IT5|). To obt ain numerical 
estimates for these types of fluctuations it is neces- 
sary to have the value of a for thin layer of TCL2O5. 
In the existing publications we have found very wide 
range: from a Ta2 o 5 - -(4.43 ± 0.05) x 10~ 5 K _1 [T3j 
(the scheme of measurement is described in ^31) to 
ft-Ta 2 5 — 3.6 x 10~ 6 K _1 (without measurement error, 
unfortunately) [Ej and a Ta2 o 5 - (5 ± 2) x 10~ 6 K" 1 
16 . It can be explained not only by possible errors 
of experiment but also by the fact that properties of 
tantalum pentoxide may be strongly depends on pro- 
cedure of layer deposition on substrate • So we use 
for estimate the value aTa 2 o 5 — 5 x 10~ 6 K _1 because 
coating measured in paper was deposited on thin 
fused silica plates by the same technology as coating 
on LIGO mirrors. The other material parameters are 
listed in Appendix IB1 

It is also useful to present estimates for thermore- 
fractive noise |Hj: TD temperature fluctuations pro- 
duce fluctuations of coating layers refraction indices 
ni and rij due to its dependence on temperature: 



Pi 



dT T ' 



P2 



^0, 
dT T ' 



(21) 



Refraction indices fluctuations produce in turn fluctua- 
tions of phase of the wave reflecting from mirror. This 
noise can be recalculated into equivalent noise displace- 
ment of mirror 1 : 



-TD 
'trefr 



P 



V2|3 2 A 2 kT 2 
7tr 2 l v / pT T Kcu ' 



1X11X2(13, + (3 2 ) 



(22) 



(23) 



4(nf-n5, 

For estimates we use "optimistic" (i.e. smallest) value 
of (3 Ta2 o 5 - 2.3 x 10~ 6 ET 1 18 (for example, [TJj 
reports unexpectedly high value |3Ta 2 o — 1.4 x 10~ 4 
K-i). 

We also want to estimate the noise associated with 
the mirror material losses described in the model of 



ifiere we present formula correcting an error in formula (3) 

m 
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Figure 2: Dependence of different noises on frequency 
in gravitational wave antenna with fused silica mirror 
and Ta205+Si02 coating presented in terms of dimen- 
sionless metric Notations: "SQL" — Hsql(cu) 

(|26fl . "Brown" — Brownian fluctuations caused by 
structure losses and described by formula Q24J1. "TD 
layer" - surface fluctuations caused by TD temperature 
fluctuations (thermoelastic) l|14l— I18|l in coating, "TD 
bulk" — the same fluctuations in the mirror volume 
@, "TD refr" — thermorefractive noise (non- 
structural damping |19j . below we denote it as Brow- 
nian motion of the surface. In this model the angle 
of losses cf) does not depend on frequency and the fol- 
lowing formula is valid for its spectral density in the 
infinite medium PUHTIH^]: 



4k R T(1 -V) 



w 



w 



27TEr 



4>. 



(24) 



Note that we have used the new value of cf) ~ 0.5 x 1 _ 
for the loss angle in very pure fused silica |24j . 

The gravitational wave antenna spectral sensitivity 
to the perturbation of metric h(tu) may be recalculated 
from the displacement X noise spectral density using 
the following formula: 



h(w) 



L 



(25) 



where we use No = 4 due to the fact that antenna 
has two arms (with length L) with two mirrors in each 
one. However, for TD layer fluctuations we use No — 2 
because in LIGO interferometer arms the Fabry-Perot 
resonators end test mass (ETM) coating TD fluctu- 
ations are considerably larger than ones in the input 
test mass (ITM) coating due to the coating thickness 
in ETM is approximately ~ 6 times larger j21) . 

The LIGO-II antenna will approach the level of 
SQL so we also compare the noise limiting sensitivity 
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with this limit in spectral form 
h-SQL = 



8ft 



mcu 



2 L 2- 



(26) 



For estimate we use the set of parameters for fused sil- 
ica mirror with TQ2O5 + S1.O2 coating using the param- 
eters listed in Appendix El The results are presented 
in fig. El 

We can see that in fused silica mirror the TD bulk 
fluctuations are negligibly small but TD fluctuations 
in layer of intcrferomctric coating are about Brownian 
one and only 2 times smaller than standard quantum 
limit at cu ~ 2n x 100 s _1 . Note that thermoelastic 
fluctuations in layer weakly depends on frequency ~ 
tu -1 / 4 (see dJ) and for frequencies > In x 300 s _1 
these fluctuations dominate completely. 

For comparison we present in fig.|31the estimates for 
sapphire mirror with the same Tc^Os + Si02 coating. 
We can see that for sapphire mirror the contribution of 
TD noise in layer is smaller and only at tu > In x 10 
s _1 it is close to SQL. 

We can see that TD fluctuations in the same coating 
Ta205 + Si02 on sapphire and on fused silica differs ap- 
proximately by 3 times. It can be explaied by the fact 
that factor k s Cs (denominator in H14|0 for sapphire is 
about two orders larger than for fused silica. Note that 
equivalent thermal expansion coefficients described by 
formulas (|16I -I18 |I are practically the same 



a, 



Ta 2 05 + Si0 2 
Si0 2 

v Ta 2 5 + Si02 
*A1 2 3 



2.8 x 10~ 6 KT 1 , 
-2.8 x 10~ 6 KT\ 



It is also worth noting that surface fluctuations 
caused by TD temperature fluctuations (thermoelas- 
tic and thermorefractive) in layer (coating) depends 



on laser spot radius weaker (~ 



r ) than " bulk" noise 



5 Conclusion 

The slow dependence of interferometric layer (coating) 
TD noise on frequency cu _1 ^ 4 and on the radius of the 
beam spot Tq 1 is worth noting. 

We can see that TD noise in interferometric layer 
is close to the standard quantum limit and it may be 
an obstacle for interferometric gravitational antennae 
(projects LIGO-II and especially LIGO-III). 

The resume of this paper may be formulated in the 
following way: it is important to measure in situ the 
value of effective thermal expansion coefficient a for 
interferometric multilayer coatings of high quality mir- 
rors for presision measurements. 
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Figure 3: Dependence of different noises on frequency 
in gravitational wave antenna with sapphire mirror and 
TCI2O5 + SiC>2 coating in terms of dimensionless metric 
Notations are the same as in figEl 
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A TD fluctuations of surface of 
half space. Langevin ap- 
proach 

Here we present the derivation of formula ljl4|l using 
Langevin approach. We calculate the surface fluctua- 
tions (averaged over the spot with radius To) of half- 
space covered by layer with thickness d. The surface 
fluctuations are originated by TD temperature fluctua- 
tions which in turn produces deformations due to non- 
zero thermal expantion. We consider the case when 
thermal expantion coefficient a is non-zero only inside 
layer and other thermal and elastic parameters are the 
same in layer and half-space. We also assume that fol- 
lowing strong unequalities fulfills: 



d<r T <r , r T 



pCcu 
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A.l TD temperature fluctuations 

We use Langevin approach proposed by yQ using ther- 
moconductivity equation for half-space z > with 
thcrmoisolated surface with fiuctuational forces F in 
right part: 



One can find solution as sum (see for example sec. 3.4 
in |H]): v = v' a ' + Vcp. Here function cp fulfills Poison 
equation (it is derived from lj2U): 

1 + y 



Acp 



1 --v 



-a du6(z — el 



(35) 



— a 2 Au = F(t, t), a 2 = — — , (27) without boundary condition and function v' a ' fulfills 



3t 



k B T 2 K 
(PC) 2 



PC 

8(t-ti) x (28) 



to H31f) with zero right part and the following boundary 
conditions: 



xAf5(x-xi)6(y--y 1 ) x [6(z - z- l ) + 6(z + zi )] ) . 

Here we introduce mirror fiuctuational forces (last term 
in square brackets) in order to reformulate problem for 
the whole space. 

After space and time Fourier transform we can 
write down the solution 



1 +-v 



9 z <p 9 2 cp~ l 



dy 2 9x 2 
9 2 cp 



(36) 



fl +"v) 9x9z 
E 9 2 cp 



'yz. 



(1 +-v) 9y9z 



u(f,t) 

u(co, k) 
(F(£,u>)P(ici,tui)) 

F 2 T 2 



dk duo 



-oc W 

F(k,cu) 
a 2 (k) 2 +icu' 



u(u>, k)e 



ia)t+ikr) 



(29) 



By this way for function v' a ' we obtain the problem of 
half space deformation with stresses (13611 on boundary. 

We are interesting of surface displacement along z 
axis averaged over spot with radius To = V2/uo in 
center of coordinate system: 



= 2 



F^T^ (27t) 4 |k| 2 6(cu - cut) x 
x5(k x - k x i ) 5(k y -k y i) x 
x[5(k z -k zl ) + 6(k z + k z ,)] 
k B T 2 K 



Xft] 



1 



dy dz 



9cp 
"9z 



x(37) 



x exp 



(PC) 5 



(30) 



x 2 +y 2 



A. 2 Equation of elasticity 



Solution of Poison equation (|35|l is known 
a(1 +-v) 



The problem of elastic deformations V for half-space 
covered by thin layer caused by thermal expansion due 
to TD temperature fluctuations u can be described by 
equation [231 125j with zero stresses cx zz , cr yz , cr xz on 
free surface z = 0: 



<Ptx,y,zj 



471(1 -<v) . 



dx'dy' x 

X) 

2u(x>,y',z<) 



(38) 



^/(x-x') 2 + {y-y') 2 + : 



1 --v 
1 +t 



x grad div v - 
1 -2-v 



2(1 +y 



■rot rot v = aVu, 



v 



1 - 2-v 
— au + 
E 



9v x 3v, 



9v, 



1 -f-v 
1 — 2-v 9v 



H + 

9x dy 9z 



Here we add solution by mirror source (~ 5(x+e)) — it 
produces the multiplier 2 in last fraction. Such function 
(31) cp makes zero contribution into surface displacement 
(see (|37H . As addition the tangent stresses cr xz and 
<jy Z in boundary conditions (|36(l are also become zero. 
Only normal pressure is non-zero and is equal to: 



adE(9 



'yz 



1 +v 9z 
2(1 +-v) \ dy 9z 



9v x 



2(1 +-v) V 9z 



9v z 
"9x 



z=0 



z=0 



z=0 



= 0, 
0. 



(32) 
(33) 



2tt(1 -y) 



dk x dk v dk z dcu 



(27t) 4 



^ F(k)T(a>) ^mt+ii^Tj+jicx x 



Q 2 |k| 2 +1CU 



Using condition d <ti To we can present the right part 
of l|31|l as following 



dy'dx' x 

— oo 

gik y (y-y ')+ik x 



aVu =^ adV(u6(z-e)) 



(34) 



v /(y-y') 2 + U-x') i 
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We take the last integral (over dy 'dx') going into cylin- 



dric coordinate system with notation Kl = y ky + k 2 
(using formulas 2.5.24.1 in |2E] and 2.12.4.28 in |S7j): 



Ii (kl_) 



p ik y («— -y ')+ik x (x'-x) 

dy'dx' 



In 



v/(y-y') 2 + (x-x') 2 

•oo gik^rcoscti £7t 

rdr = — — — 

o vr 2 



Gathering we obtain expression for cr 
ocdH 



(1 -y) 



dk x dk y dk z dcu 



x k ± F(1c)T(a>) eicut+ilCyU+ilCxX (39) 
a 2 |k| 2 +ia> 

Now we can write down the expression for v' a ' at z = 
(see sec. 8 in [2"5]h 



(a) _ 



1 -V 



dy'dx' 



a*=°(y',x') 



v/(y-y') 2 + (x-x') 2 

Substituting it into l|37|l we obtain: 

dk x dk v dk z dcu 



X(t) 



2ad(1 +^ 
F(k)T(cu 



(27t) 



(40) 



(a 2 |k| 2 +ia>) 

One can write down the time correlation function (us- 
ing 



BfT) = (Xft)X(t + T)) = 



2 d 2 (l +^) 2 


>00 






— OO 


>oo 


r0 ° dk z k ± dk 


J — oo • 


o (2tt) 2 



27T 



1 '0 1 



(41) 



Now we can write down one side spectral density: 
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Sx(cu) = — (1 +^) 2 a 2 d 2 F 2 T 2 x 

,2 , V 2 



■°° (k 2 +k 2 )dk z k ± dk ± _ k 2 r 
(Q 4 (k 2 +k 2 ) 2 + 0» 2 ) 6 



(42) 
^/ 2 
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b 4 = 



2r a 4 J 

2 2 



ly + y/y 2 +b 4 
(y 2 +b 4 ) 



e-« dy, (44) 



4a 4 



B FDT approach calculations of 
thermo-elastic noise 
in interferometric coating 

Here we derive formula Ijl4(l using fluctuation- 
dissipation theorem (FDT) 28 29 by the following 
thought experiment: 

We imagine applying a sinusoidally oscillating pres- 
sure, 



-r 2 /ri 



7trt 



(45) 



to one face half-infinite mass covered by layer. Here 
F D is a constant force amplitude, u) is the angular fre- 
quency at which one wants to know the spectral density 
of thermal noise, and the pressure distribution i|45|) has 
precisely the same spatial profile as that of the gener- 
alized coordinate X, whose thermal noise S x (tu) one 
wishes to compute. 

The oscillating pressure P feeds energy into the test 
mass, where it gets dissipated by thermoelastic heat 
flow. Computing the rate of this energy dissipation, 
Wdi S s, averaged over the period 2n/a> of the pressure 
oscillations we can just write down (in according with 
fluctuation-dissipation theorem) the spectral density of 
the noise Sx(tu) is given by 



SxM = 



8k B TW, 



diss 



F 2 < 



(46) 



The rate Wdi ss of thermoelastic dissipation is given 
by the following standard expression (first term of Eq. 
(35.1) in H3): 



T 



(Vu) 2 rdc|) drdz 



(47) 



Using condition a /Tq <C w (and assuming that k ~ 
k 2 ) we can easy take integral over kj_ then over k z and 
obtain the formula l|14|). 

Note that formula l|42|l can be used for calculation 
Sx with arbitrary value cu: 



S x (u>) = ^(1 +^)Vd 2 F 2 T 2 x I X 



(43) 



V2r a 4 



(x 2 +y) dxdy 
(x 2 +y) 2 +b 4 ) 



Here the integral is over the entire test-mass interior us- 
ing cylindric coordinates; T is the unperturbed temper- 
ature of the test-mass material and u is the tempera- 
ture perturbation produced by the oscillating pressure, 
k is the material's coefficient of thermal conductivity, 
and (. . .) denotes an average over the pressure's oscil- 
lation period 2tc/ w (in practice it gives just a simple 
factor ((me iwt ) 2 ) = 1/2). 

The computation of the oscillating temperature 
perturbation is made fairly simple by two well-justified 
approximations 1 : 

First: Quasistatic approximation. We can approx- 
imate the oscillations of stress and strain in the test 
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mass, induced by the oscillating pressure P, as qua- 
sistatic. 

Second: We assume that the following strong un- 
equality is valid 



2(1 +T b )(1 -2v h )P 



3 z u z 



(1 +v b )(l -2v h )P 



(52) 
(53) 



d <C r T <C r D , r T = 



pCcu ' 



here tj is length of diffusive heat transfer and p is the 
density and C is the specific heat . This means that, 
when computing the oscillating temperature distribu- 
tion, we can approximate the oscillations of tempera- 
ture as adiabatic in transversal direction (tangent to 
surface) but non-addiabatic in normal to surface direc- 
tion. 

The quasistatic approximation permits us, at any 
moment of time t, to compute the test mass's internal 
displacement field v, and most importantly its expan- 
sion 

6 = div v, (48) 

from the equations of static stress balance (Eq. (7.4) 
in gSj) 



(1 -2t)V z v + V(V • v) 







(49) 



(where "V is the Poisson ratio) , with the boundary con- 
dition that the normal pressure on the test-mass face 
be P(r,t) [Eq. l|45|l] and that all other non-tangential 
stresses vanish at the test-mass surface. 

Ones has been computed, we can evaluate the 
temperature perturbation u from the thermoconduc- 
tivity equation 



(9 t -a 2 A)u = 



-aET 



pC(1 -2v 



■3 t e, a 2 = — (50) 
pC 



Here E is Young's modulus. It worth to underline that 
we can not use addiabatic approach in l|5U|) due to small 
thickness of layer. 

This temperature perturbation u can then be 
plugged into Eq. H47|) to obtain the dissipation Wdi SS 
as an integral over the gradient of the expansion. This 
Wdiss can be inserted into Eq. Q46(l to obtain the spec- 
tral density of thermoelastic noise. 

B.l Elastic problem 

Elastic infinite half space. We can assume that 
layer does not influece on deformations in substrate due 
to its small thickness. Then we can use the solution to 
the quasistatic stress-balance equation (|4^|) given by 
a Green's-function expression (see (8.18) in (2H|) with 
Fx = =0, F z = P(r), integrated over the surface of 
the test mass: Below we will need the expression for 
0' b ' = divu. Using results of PP one can calculate 
longitudial and transversal parts of 0' b ' separately: 



e[ b) = 



(3xUx + dyUy) z = 



(51) 



Layer. For solution of thermal problem we will need 
to calculate expansion O' 1 ' in layer. We assume that 
deformations of layer in transversal plane are the same 
as in substrate (the same natural assumption was done 



in i- e - ©II = 0[| U 'lz=o- One can use equation 

(5.13) for stress in 23 for calculation @± = it 



(b). 



-P = 



Ei 



(1 -2vi 



(54) 



-<vx (u( f '+u' f '; 



Using this equation one can find u zz and full expansion 
0i in layer: 



u 



i) _ 



©i 
Y b 



Yi(1 -vi) 
P 

"Yi(1 - -vi) 
E b 



1 - 



(1 +-VbH1 -2*i 



•viYi 
Y b 

Yi(1 - 2vtj 
Y b 

, Y x = 



(55) 
(56) 



1 -2vi 



B.2 Thermal problem 

The expansion in layer and in substrate are the source 
of heat with power per unit volume (see for example 
sec. 31 in [23]) which are equal to Wi inside layer and 
to Wb in substrate 



W l = ? =^Iia J l( 
(1 -2*i) 



-a b E b T . 
1 -2-Vb) 



It seems that we can divide our problem into two prob- 
lems: 

First: The "bulk" problem. In half space there is 
power source Wb . It corresponds to the problem solved 
mill. 

Second: The "surface" problem. There is power 
source only inside layer. Below we consider only second 
problem. 

We have thermal conductivity equation in layer 
with constant right part. Using conditions (jSJ one can 
simplify thermal conductivity equations to one dimen- 
sion form as following 



1 



Kb 



1 - 



itu Cb 
icu Ci 



di 6Tb = 0, 



3: 6Ti 



(57) 
(58) 
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w 



A 



(l-2vi)d (1-2v b )C b 
2(1 |7 b )Ta b 



PA, 

1 +"Vl 



C b 

oci C b 

a b Ci 2(1 — -vi)(1 +^/ b 
, Ei(l -2v b )(1 +^ b 



(59) 
(60) 



E b (1 +v t ) 

We have condition of thermal isolation at z = and 
condition of continuity of temperature and heat flux at 
z — d. So we can find the solution of system (|57l I58fl 
as following 

1 +i 



STi = A — Bcoshyiz, Ji = 



Kl 



6T b = Ce- yb(z - d) , 



Tb 



V2 

1 + i / w C b 



\fl V K b 

C = A-Bcoshyvd, (61) 
K bYbC = Ki Yx B sinhyvd, (62) 

wR sinhy^d 



A 
B 
R 



w, C 



R sinhyvd + coshyid' 
w 



R sinhyid + coshy^d ' 

KlYl 



(63) 
(64) 
(65) 



K b y b 

Due to conditions JJjJl and hence |y g d| <C 1 we can 
simplify espressions for temperature: 

6T\ ~ w(1 — cosh yiz+ Ryid cosh Yi.z), (66) 
w KiY\ d 



6T b 



-Tb (z-d) 



KbYb 



(67) 



We are interesting only in temperature distribution in 
infinite space because volume of layer produces small 
contribution into total budget of dissipated energy: 



W diss = 



c 2 

cu 2 d 2 T 



(68) 



^VCbKbtU 

Using this result and 146(1 one can write down spectral 
density 

4v/2(1 +^ b ) 2 k B T 2 2 2 

oc off d , (69) 



S x (tu 



7rr 2 VC b K b w 

Oeff = «b A, 



(70) 



B.3 Generalization for multilayer coat- 
ing 

We can generalize this result for multilayer coat- 
ing consisting of N alternating sequences of quarter- 
wavelength dielectric layers. Every odd layer has the 



refraction index ni , thickness di = A/4n.i (A is the 
incident light wavelength), its parameters we mark 
by subscript i . Every even layer has parameters ri2 , 
&2 — A/4ri2, its parameters we mark by subscript 2- 
It can be easy shown that the equivalent thermal ex- 
pansion coefficient a c g of total coating and the total 
thickness d of coating are the following: 



Ocff 

d 
A 1 



A 2 



a s di a s d 2 . 

Ai + - — — — A2 , 



di + d2 di + d2 

N(di +d 2 ), 
Cfi QCfi x 
C b 2a b 

1 +-Vfi 
1 --v f i)(1 +^ b ) " 

Cf2 . <Xf2 
1 X 

Cb 2ab 

1 +^f2 



E fl (1 -2-v t 



Eb(i -Vfi ; 



(71) 
(72) 
(73) 



(74) 



-f2 



1 -2v* 



(1 --Vf 2 )(1 +"Vb) Eb(1--V f2 ) 



And the spectral density Sx(cu) is defined by substi- 
tution (|71|1 into formula (|69|) . We have checked it by 
solution of thermal problem and calculation W<ji ss for 
one pair of layers, for two and three pairs of layers. 



C Finite sized mirror 

Our calculations for finite sized mirror are based on re- 
sults of Liu and Thorne J2| ■ The order of calculation 
steps is the following. First step: we calculate total 

f s I 

expansion S and its tranversal part in substrate. 

Second step: we substitute calculated value 01, into 

(|54fl and calculate Q± and then 0f. Third step: we 
substitute the values O s and Of into thermal cunduc- 
tivity equations (|57l and lSSj l. Forth step: we calculate 
dissipated energy Wd; ss - 

Calculations of O s . We use the results and nota- 
tions of ^2]- The radius of mirror is R, its height is 
H. 

8 S = Foe** f_ P°n-2. s ) + 2(1 -av.)co ^ 5) 



2jk m A m + B m ] J (k m r) J , 
m / 



A = 



Ev 



(1 -2v)(1 +y] 



2(1 +V)' 



(76) 



1 r 6R 2 f Jo(C m )p m (77 , 

m=1 111 

Here Cm are the m-th root of equation Ji (x) = 0, 
JoM> JiM are Bessel functions of zero and first or- 
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ders. 

Pm 

B m 
D m 

Qm 

Pm 
Ym 
6m 



Pm 



2(1 -T s )c 



exp(-k m r 2 /4) _ Cm 

7tR 2 J 2 (Cm) ' m R ' 

A m | z= o = Ym + 5 m , 



= d z B m (z)| z=0 = k m • 



(Tm + § 
[km A m + B' ] = 



|X s (6 m - ^ 



exp(-2k m H) 



Pm(A s + 2u. s ) 


1 - Qm 


+ 2k m HQ m 


kmm(A s + m) 


(1 - Qm) 


2 -4k 2 n H 2 Q m 


p m (A s + 2u s )Q 


1 - 


Q m +2k m H 


k m Ll s (A s + H s 


(1 - Q 


m )2-4k2 a H2Q m 


-Pm 




+ ^s(1 - Qm) 


2k m Li s (A s + p. s ) 


(1 - Qm 


) 2 -4k m H 2 Q m 


-PmQm 


1m ~ 


Ml - Qm) 


2k m Li s (A s + lx s ) 


(1 - Qm 


)2-4k2 n H 2 Q m ' 



D, 



2k m H 2 (A s + H s )±2p. s k m H, 
Pm (1-Qm) 2 



(A s + ^ s ) (1 -Q m ) 2 -4k m H 2 Q m 



Calculations of 



e 



(s) 



Foe icot ('2-v s po , 2(1 -v s )c 



+ (78) 
(79) 



+ Y_ k m A m |z=o Jo (k m r) J , 

m / 
km(S m + Ym) = F X (80) 



(1 -Qm) 2 d -2T s )+4Q m k m H 2 
(1 -Qm) 2 -4k m H 2 Q m 



Calculations of expansion 0f in layer. Using 
we obtain: 



u z 



-P(1+^ f )0-2v f ) . e (.)_V!_ |(81) 



E f (1 --v f ) 



II 1 



Tf 



-P(1+-vf)(1-2v f ) +ef . ) l-2v i 



E f (1 --v f ) 



II 1 



"Vf 



It is useful to present pressure P as series and rewrite 
Of as following; 



P = ¥ e iwt (po+ X PmIo(k m r 
Of 



(83) 



m=1 



< - 2v f iwt . 

to e po 



1 -"Vf 



2t s (TjPVfj 
E s E f 



+ 



„ (1+Vf) , (1+-v.. 
K m = — - H = x 



-XPm!o(k m T)K m ,(84) 

m=1 / 

ri 4-1/ i 

(85) 



E f E_ 
(1 -Qm) 2 d -2v s )+4Q m k m H 2 



(1 -Q m ) 2 -4k 2 ,H 2 Q m 



Calculations of right part in thermal conductiv- 
ity equation. Using l|58|l we obtain the expresion for 
w: 



w 



-a f EfTO f a s E s T0 s 



CfW 



(l-2v f )Cf (l-2v.)C,' 
a s TE e ia,t x 



(86) 
(87) 



C f a f 
-Pohr + - 
C s a s 



2-v s E f 1 +-v f 



(1 --Vf)E s 1 --Vf 



•Cf a f E f (1-T 



C s a s E s (l --v f ' 

OO ^ 

Y_ PmJo(kmr)L 



L m = 



m=1 

«f (1 +-v f ) 1 + t s 



a-s(1 -Vf) 1 — -Vf 

(1 +VsK1 -Qm) 2 

(1 -Qm) 2 -4k 2 R H 2 Q m ' 
'a f E f (1-2v s ) 1 

X [ + 



2Cf 



a s E s (l -T f ) 1 -^f C s 



(88) 



Calculations of dissipated energy Wdi SS . Using 
(|4"7jl we obtain: 



v = 



a 2 F 2 Ta) 2 d 2 
5_JJ -p 

2\/2 7tR 2 v / k7C 

ocf f\ +"Vf 2-v s E f 



oc s V 1 -^f n -v f )E s 
C f afE f (1 --v 



C s a s E s (l -T f ) 



" c s 

n 2 



12R 2 

"H 2 "^ C 2 J (C r 



m=1 



s 2 = £ 



-<i-o/ 2 



m , J§(Cr 



(89) 

S 2> 
(90) 

(91) 



Here we used the (nonstandard) orthogonality rela- 
tions: 



Jo(kmr) J (lCnT)rdr 



R 2 
2 



Jo (Cm) 6 m n , 



J (k m r)rdr = 0. 
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Now one can calculate numerically the coefficient References 



c 



ism 



W fsm 
diss 



2R 2 (1 +^ S ) 2 A 2 



(92) 



using flH2) and 

We numerically estimate Cft m for R = 19.4 cm, H = 
1 1 .5 cm, To = 6 cm and for TCI2O5 + S1.O2 coating on 
fused silica substrate: 



CftE 



D Parameters 



m 

a 
P 
E 



dE 
dT 



E 

dE 
dT 
n 



E56 



= 27tx100s~\ T = 300 K, r ~ 6 cm 
= 3x10 4 g, A = 1.06 u, E = 4x10 5 cm;; 



n = 



a 



a 
a 
a 
E 



C = 

n 
tin 

dT 
dn 

dT 
dn 

dT 



Fused silica: 

y-7 



5.5 x ICr' K~ , k = 1.4x 10 

2.2 -§L C = 6.7x 10 6 -^1, 
cm- 4 g K 

7.2 x lO 11 ^, v = 0.17, 

cm 3 

0.5 x 10~ 8 El; 

8 er S 



5 er § 



cm s K ' 



— = -1.5x10 



K cm' 



3 ■ 



I.45, — = 1.5 • 10~ 5 K _1 , 
dT 

Sapphire: 

5.0x10- 6 KT 1 , k = 4.0x10 6 erg T . 

cm s K 

4.0 -JL C = 7.9x 10 6 

cm- 4 g K 

4x10 12 -^4, v = Q29, ct) = 3x10- 9 , 



_ = -4x10 



cm J 
8 erg 
K cm 3 



1 .76, (A = 1 urn 1521 
Tantal pentoxide TCI5O5 
-4.4 • 10~ 5 KT 1 HI, 

3.6 • 10 6 kt 1 nni, 

(5± 1) • 10~ 5 K" 1 M 

1.4 x 10 12 -^ 30, ^ = 0.23, H3 

cm 3 

4.6 x 10 6 ° rS . , p = 6.85-^, 

cm sK cm 3 

3.06 x 10 s 2£ 

g K 

2.1, 

1.21 ■ 10~ 4 K _1 (A ~ 1.4am) ED 
2.3 • 10~ 6 K- 1 (A ~ 0.63am)HB], 
4.7- 10~ 6 K _1 ED 
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